Quark Schwinger-Dyson Evaluation of the 11, 12 Coefficients in the Chiral Lagrangian. 
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Using a systematic expansion of the quark-antiquark Bethe-Salpeter wavefunctions in the rela- 
tivistic quark model and working to 0(P 4 ), in the chiral limit, we are able to derive theoretical 
expressions relating the coefficients of the chiral lagrangian l\, h to the underlying quark-antiquark 
wavefunctions and interaction kernels. This is accomplished by using a novel technique based on 
a Ward Identity for the quark-antiquark ladder kernel which greatly simplifies the required effort. 
Numerical evaluations are performed in two simple specific models. 



I. INTRODUCTION. 

It has traditionally been considered a triumph of the- 
oretical physics when the parameters of an effective, low 
energy theory which correctly describe phenomena at 
a given scale can be related to those of an underlying, 
more fundamental scheme of thought which grounds it. 
Brilliant examples are Fritz London's explanation of the 
quantum nature of the Van der Waals forces [l| or the 
derivation of the atomic relativistic corrections as a con- 
sequence of the Dirac's equation for the electron. Low 
energy hadronic processes are interpreted with the aid of 
mainly two types of theories: nucleon-nucleon non rela- 
tivistic interactions such as the Nijmegen Q or Argonne 
potentials, for the heavier hadrons, and relativistic 
chiral lagrangians Q| for the lightest components, the pi- 
ons. 

The deeper quark theories such as QCD or any mi- 
croscopic models thereof pretend in principle to describe 
the totality of hadronic physics. They attempt to be 
complete descriptions of hadronic processes. Unfortu- 
nately, the complexity of many body hadronic calcula- 
tions makes it forbidding to fully exploit the underlying 
scheme, and maintain the validity of the low energy ef- 
fective theory. 

As a consequence, an initial goal for the microscopic 
theory should be to reproduce in some limit the macro- 
scopic models and to relate their parameters to its own 
set (hopefully smaller). In this paper we make the case 
for microscopic quark models inspired in QCD as gen- 
erating the parameters of the chiral lagrangian. This 
lagrangian, describing the low energy behavior of a pion 
system, and being able to incorporate the coupling of pi- 
ons to other mesons (as much as the low energy theorems 
ofPCAC [| do) is universal (in the sense that any theory 
with the same symmetries can be cast in its form) and 
provides a consistent derivative expansion in powers of 
the momentum and mass of any pions present in a sys- 
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tern, divided by a typical scale of the strong interactions. 

Unfortunately, this derivative expansion has to incor- 
porate new coefficients order by order. These new co- 
efficients absorb the divergences of loops generated by 
the vertices of smaller order terms, and so their value 
is generally renormalizcd. Still, the common usage of 
this Lagrangian ;6j proceeds by fitting this coefficients to 
some observable set at a given scale. 

We show how these coefficients can be related system- 
atically to quark level parameters in the planar approxi- 
mation. This has been accomplished in the past for the 
simplest, 0{P 2 ) chiral lagrangian whose parameters are 
only two, in usual notation, M n , / w (the pion mass and 
decay constant). To this order, these two parameters are 
conventionally set to take their physical value. To the 
next order, the lagrangian contains six parameters, which 
generate the 0(P 4 ) vertices, h, I2, absorbing divergences 
in the 4 pion Green function, £3, 14 which absorb countert- 
erms of the mass and axial current renormalizations, and 
finally M ff , f^. The complete renormalization scheme is 
specified in |4j. The parameters M ff , / w have long been 
accounted for by relativistic quark models The Vs 

on the other hand have not been treated in quark mod- 
els with non-contact interactions. We insist in the point 
that any theory which respects the SU(2)l x SU(2)r 
chiral symmetry breaking pattern, let it be a Nambu- 
Jona-Lasinio quark theory [!j , a large N c expansion |10| , 
a string theory or any other exotic creation, can be cast 
in the form of the chiral lagrangian and the only differ- 
ence between all of them is the numerical values of the U 
coefficients. 

It is therefore of paramount importance to determine 
them from the theories which we believe correctly de- 
scribe the physics at the GeV scale, in terms of quarks 
and antiquarks. Lattice determinations are making 
progress in that direction ^f|, but the Schwinger-Dyson 
equation formalism should provide an alternative deter- 
mination in the near future. An interesting paper |T^ 
exists where, at the lagrangian level, the action for a rel- 
ativistic quark model is bosonized to obtain an effective 
meson lagrangian, used then to calculate pion-pion scat- 
tering lengths. We are going to theoretically extend this 
approach in two directions. First, we will start with the 
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most general chirally symmetric quark model, in which 
the pion is well described by a quark-antiquark pair af- 
ter chiral symmetry breaking (encompassing in this way 
an ample spectrum of models) and by using their chiral 
properties, reduce the four pion Green's functions to a 
minimal set of diagrams. In this way, no bosonization is 
performed, and at all steps the way quark interactions 
arrange themselves to comply with the chiral theorems is 
explicitly visible. Second, comparing the result with the 
same calculation in a chiral lagrangian formalism, one 
can immediately read off the Zj coefficients of the chiral 
lagrangian in terms of diagrams which can numerically 
be calculated in the quark model. This rather technical 
numerical evaluation will be simplified in this work by 
confining ourselves to simple, finite models, although the 
numerical results will then be limited. The method used 
here has already successfully being exploited to demon- 
strate how this class of models comply with the Weinberg 
theorem in [la. Ill fl5| . The Weinberg theorem was de- 
rived with an expansion to 0(P 2 ), 0{M V 2 ). We now con- 
centrate on the <3(P 4 ), 0(M W °) chiral lagrangian, that is, 
the only parameters are f w , l±, and We will perform 
the same expansion in the quark-antiquark diagrams and 
compare the results to read off l\, h- The expansion will 
be carried out whenever possible in a Feynman diagram 
language to avoid lengthy expressions for the sake of read- 
ability. The rest of this paper is organized as follows: in 
section|n]we briefly settle the notation for our chiral per- 
turbation theory discussion and remind the reader of a 
few well-known facts in this field. Section ITTT1 settles the 
notation of the microscopic quark manipulations to fol- 
low and provides the reader with a useful chiral Ward 
Identity recently introduced 0, 01- Section [Tvl is the 
core of the paper and presents the reduction of the pion 
scattering amplitude, whereas the resulting diagrams are 
calculated in two simple models in section [V] Some is- 
sues clarifying the normalization of the Bethe-Salpeter 
equation are relegated to the appendix. 



where the scalar product dots are in isospin space. This 
lagrangian is on shell, for massless pions (else the I3, I4 
counterterms should also be present) and contributes at 
tree level to the 0(P 4 ) pion-pion scattering amplitude, 
and it is this contribution which we aim to reproduce 
microscopically. In the chiral formalism, there are also 
one-loop contributions from the 0(P 2 ) lagrangian which 
we do not consider in this work, since our quark-level 
calculation will not be extended to meson loops. There- 
fore, to this level, it is fair to compare our results only 
with those obtained in chiral perturbation theory with- 
out meson loops. With this caveat in mind, the pion-pion 
scattering amplitudes are straightforwardly obtained. By 
using crossing symmetry, the different isospin channels 
can be related in terms of only one amplitude A: 

Ti =2 =A(t,s,u)+A(u,t,s) (3) 
T1—1 — A(t, s, u) — A(u, t, s) 
T/ =0 = 3A(s, t, u) + A(t, s, u) + A(u, t, s) . 

This amplitude A(s,t,u) can be obtained from the pro- 
cess 7r+ 7r_ — > 7To7To. Due to the final state Bose symme- 
try, and in the chiral limit when the Mandelstam vari- 
ables satisfy s + 1 + u = 0, the most general amplitude of 
order P 4 containing the polynomials s 2 , t 2 , u 2 , st, su, tu, 
reduces to Ais 2 + A^ii — u) 2 . The coefficients obtained 
from the lagrangian above yield 

A^(s, t,u) = jr 4 [(2h + l j)s 2 + l j(t u) 2 } . (4) 

A full discussion of this and related issues, for example 
the relation between f n and F which we further ignore 
in this paper to the order we are working can be found 



III. NOTATION FOR QUARK MODELS AND 
CHIRAL WARD IDENTITIES. 



II. CHIRAL LAGRANGIAN OF ORDER P 4 . 

The macroscopic theory one generally writes down for 
pion fields alone is to lowest order the Non Linear Sigma 
Model. One can proceed by constructing, from the three 
pion fields, 7? — (tv\, 7r2, ^3), a 4- vector normalized to one 
(this normalization is equivalent to eliminating the ex- 
plicit a degree of freedom from the Linear Sigma Model) 
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and then constructing Lorentz scalar, parity invariant 
terms. To 0(P 4 ) that la grangian can be extended by 
terms which in the chiral limit (m q = 0) have to be of 
the form Q 

1 



A pion with momentum P couples in relativistic mod- 
els to fermion lines whose momenta will be denoted by 
k, k' . In the massless quark limit, whenever P = 0, then 
k = k' . We start by considering the bare fermion propa- 
gator from any standard quark theory, 



S (k) = 



ft — m + ie 



(5) 



and, after spontaneous chiral symmetry breaking medi- 
ated by a strong interaction 



the full fermion 



(1) propagator parameterized as 



S(k) 



A{k 2 ) /t-B(k 2 )+ie 



(G) 



which we take to be a solution of the planar Rainbow 
Schwinger-Dyson equation 



(2^ 



V a S(k + q)V a K(q) . (7) 
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We further define the bare axial vertex which couples a 
quark-antiquark pair to a pseudoscalar current by means 
of the shorthand Ja (notice that m u = = in the 
chiral limit employed in this paper): 



7l = y H^7 M 75 + 2im u75 ) 



(8) 



which satisfies 



TX(M') - T (^ 1 (fc)75 +75% 1 (fc')) (9) 

and the dressed axial vertex, dressed with a planar ladder 
T A , given by 

T A (k,k')= 7 %(k,k') + J V-Sih + q^Si^ + ^VaKiq) 

( 10 ) 

or reconstructing the planar ladder expansion (in graph- 
ical form): 



=0= 



I|1U|) takes the form 



ST A S = 



(11) 



(12) 



from which one can deduce the axial vector Ward Iden- 
tity: 



(13) 



This is analogous to the abelian vector Ward-Takahashi 
Identity which in terms of the vertex defined by 



ST a u S 



In 



yields 



i(k„ - K)T^(k, k') = s-\k') - s-\k) . 



(14) 



(15) 



Next we introduce the bound state formalism for quark- 
antiquark systems. To this end we remind the reader 
of the Bethe-Salpeter (BS) amplitude \ (see 0, E3) 
for further details) which satisfies a homogeneous Bethe- 
Salpeter equation: 



X b (P,k) 



or in graphical form: 
k-P 



V a S(k'+^) x b (P,k')S(k'-^)V a K(k-k') . 

(16) 



k-P k'-P 



(17) 



Each incoming or outgoing pion in a particular pro- 
cess must contribute with one of these \ functions, which 
carry pseudoscalar quantum numbers by construction 
|19|. The BS amplitude for a particular pion depends 
on the total momentum of the pion P, and the momen- 
tum of its fermion component k ± P/2. Notice that this 
equation is the homogeneous part of l|ll)fl above when we 
interpret the pion momentum P as k — k! in the vertex 
definition. 

Now let us deepen our study of the vertex T A . From 
(|13fl can easily be seen that in the chiral limit (m q —0), 



T a A (k,k' = k) = 2iB(k 2 ) l5 - 



(18) 



in terms of the SD amplitude B solution of (7J|. Equa- 
tions l|16(l and l|10(l . homogeneous and not homogeneous, 
coincide when "f A = 0. This is satisfied in the limit 
rn q = when also P = as can be seen explicitly from 
@ and allows us to identify, up to a normalization con- 
stant, Xtt{P — 0, k) with TA(k,k' = k). This constant 
coincides with the pion decay constant in the chiral 
limit (the proof is sketched in appendix ^) and finally 
entails, in combination with 118(1 . 



xliP^^k) 



-iT a A {k,k' = k) _ 2B(k 2 ) a a 



75- 



(19) 



In [2fJ the proof was given why this BS amplitude, in 
connection with the Axial Vector Ward Identity makes 
the pion a Goldstone boson. In terms of our notation 
this has been rewritten in p^| . 

This discussion suggests a strategy to systematically 
organize the corrections to the chiral, low momentum 
limit, in an analogous fashion to that used in Chiral Per- 
turbation Theory. Since the vertex T A and diagrams 
constructed thereof satisfy interesting Chiral identities, 
let us define 



X a (P,k) 



-iT A (P,k) + A a (P,k) 

u 



(20) 



where the function A(P, k) so introduced can be ex- 
panded in a Taylor series for low P. This expansion will 
organize the momentum corrections to any diagram. One 
will first use the chiral results for T a, which will provide 
one with exact low energy theorems, and the numerical 
corrections as P is increased can then be expressed as 
overlaps of A functions. 

We do not yet specify the color, spin, flavor or mo- 
mentum structure of the interaction kernel and vertices 
V a V a K(q), except for one property: it must be chiral 
symmetry preserving, that is, V commutes with 75. This 
guarantees the satisfaction of the following Chiral Ward 
Identity (also discussed in , Q ) which proved essen- 
tial: 



3 = 



1 • < 



s- 



7-1 



(21) 



4 



75 



757A 



or, for a general vertex not necessarily pseudoscalar: 



T A 



75 



rQr> y + rQr>^ y 



(22) 



This identity allows the reduction of terms with two 
axial vertices and is the core of the present calculation. 
(We remark that S I_1 S' = / is introduced in i|21|) and 
S completes the ladder in eq. 1)11(1 leaving the explicit 

S-K), 



A. Further Ladder Properties. 

We start by observing that the pseudoscalar ladder can 
be Laurent-expanded around its pion pole. Keeping only 
the hrst term, containing the pole, one obtains: 



(23) 



where 



p2 _ m 2 



(24) 



(in the calculations contained in this paper, = 0.) 
Combining this together with the definition of A in eq. 
l|2Tl)) one can use 



£1 



(25) 



(f-^)cc 



The ladder can also be expanded in powers of the ex- 
ternal momentum: starting from the geometrical series 
(|ll|l and expanding all propagators, then resuming when 
possible, we can show that in analogy with the matrix 
relations 



1-1/ 



= 2 



1 



1 - 



1 I „.' I 1 



i-y/ V 1 - 2// V 1 -^ 



i-i// \i - y 



one has 



/c+P 1 — 1 q+P 



ML 



-(26) 



AU4 



(fc'+-P) (2) 



ML 



(fc / +P)< 1 >(fc"+P)( 1 



where the following momentum expansion of the prop- 
agators is meant by the super indices ( 2 ) 



S{q + P/2)\ P=0 



Ajq + B 



& i S(q + P/2)\ P=0 =<y lt 



A/2 



A 2 q 2 _ B 2 



d^d^S(q + P/2)\ P=0 



A 2 q 2 _ B 2 

A' jq + B' (A jq + B)(A 2 + 2q 2 AA' -2BB') 
A 2 q 2 - B 2 {A 2 q 2 - B 2 ) 2 

jqA' A 2 + 2q 2 AA' - 2BB' 



(27) 
(28) 



+2q 2 



' A" jq + B" 
2(A 2 q 2 — B 2 ) ' (A 2 q 2 -B 2 ) 2 



""^ ' ~ '~>^-v A 2 q 2 -B 2 
A jq + B f{A 2 + 2q 2 AA' - 2BB" 2 



(AA + 2 q 2 ( J qA' + B / )l29) 



A 2 q 2 _ B 2 



(A 2 q 2 -B 2 ) 2 
- (2AA' + q 2 (A') 2 + q 2 AA" - {B'f - BB" 



(A 2 q 2 -B 2 f 



[(A' jq + B')(A 2 q 2 -B 2 )-{A jq- + B)(A 2 + 2q 2 AA' - 2BB') 



The last diagram in (|2(j|) contains an annoying explicit rung of the interaction. This can be eliminated at the 
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cost of adding a diagram to any expression where it ap- 
pears: in analogy with 



y 



l 



1-1/ 1 - y 



i 



one can use 



30r = =0= 



(30) 



IV. PION-PION SCATTERING 

A. Generalities. 

The pion scattering amplitude with the Bethe-Salpeter 
and planar approximations can be derived from the fol- 
lowing Feynman diagrams: 



X-7T1 



X^4 

C: 



X"7T2 X7T3 X7T2 



(31) 



where the first two terms provide all possible planar 
topologies, but upon substitution of their zeroth or- 
der is seen to be double counted, hence we subtract it. 
The two approximations involved are: first, coupling of 
the pion to higher Fock space states is not considered and 
second, only planar diagrams are utilized. It was shown, 
using the axial Ward Identity, that these two approxima- 
tions are consistent within the Schwinger-Dyson method 
in the rainbow approximation for the fermion mass gen- 
eration and in the ladder approximation for the bound 
state equation 0,0]- This is also consistent with past 
work 21] on resonance exchange, and is equivalent to 
the lowest order in a 1/N C expansion. Reduction of this 
combination of Feynman diagrams to 0(p 4 ), 0(M % ) is 
our goal. The calculations in this section will treat the 
P's as incoming momenta. Matching the dummy Pj to 
the incoming q^ , qi 2 and outgoing q 0l , q D2 physical pion 
momenta leads to 6 different permutations, namely 



\Pi,P2, Pz, Pa 



( fan, Qi2, ~q 2, -q a i) \ 

fail, qi2, — lol, — Q02) 
fail, —lol, —Qo2, <?i2j 
fail, ~<7o2, —<lol, <?i2j 
fail,— lol, Qi2,—qo2) 
\fail,-Qo2, %2>-?ol)/ 



(32) 



where the first momentum is fixed to avoid double count- 
ing by rotational symmetry of the tt-t: scattering ampli- 
tude J3TJ). 



We concentrate on A(s,t,u), the amplitude for 
7r + 7r_ — > tto^Oj and use the following isospin wavefunc- 
tions 



1/ 

7T+ = -((7l + ia 2 ) 

1/ 



(33) 



1 



to 



V2 



CT3 • 



From them follow the traces 

Tr(7r + 7r_7r 7ro) 



(34) 



(which multiplies the four first permutations in 132|) 
above where the two charged pions are adjacent) and 



Tr(7r + 7r 7r_7r ) = - 



(35) 



(multiplying the last two permutations in (|32J) where the 
two charged pions are in opposite corners of the ampli- 
tude). 

Finally the Mandelstam variables in the chiral limit 
satisfy 



s = 2q il q i2 = 2q ol q o2 
t = -2qnq i = -2q l2 q o2 
u = -2qnq o2 = -2q l2 q ol . 



(36) 



Since isospin has been factored out, we can ignore it 
in the rest of the calculation. Instead of working with 
r a we employ T as defined in eq. (|13fl . Accordingly, we 
define 



x := T2 X 



(37) 



to yield the normalized isospin wavefunctions in (|33|l and 
the normalization for x, Ta will be 



X 



-iT A 



V2U 



(38) 



We now start evaluating the Feynman amplitude in 
terms of the P's. Start by employing 1)38(1 to treat the 
product of four BS amplitudes: 



X7T1 X7T2 X7T3 X7T4 (39) 



\V2~U 

= (^j=T ) ^A 1 A 2 A 3 A4 + yA 2 A 3 A4 + perm. 

where the omitted terms, contain an increasing number 
of powers of T. Next we proceed to a term by term 
analysis of this expansion, further explained in the paper 

El. 
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B. Contribution with 4-As. 



Ai A 4 Ai 

ID - CZTD 

— • — — • — * — • — 

A 2 A 3 A 2 



The term with four A's is model dependent and no 
chiral properties can be used to simplify it, since it con- A4 
tains corrections to the BS pion wavefunction beyond | | + ^ * 

the zero momentum limit. Without evaluating it explic- { { As 1 — 1 

itly in a particular model yet (but see later), we can pa- 
rameterize it. To the order mP we work here, for 011- 
shell pions, if = for all i. Therefore, the 4-A dia- = 3d 1 (P 1 -P 2 Pa-Pi + Pi-Pa PrPj) + 3<2 2 P r P 3 P 2 -P 4 
grams can only be a combination of products of different 
momenta, PiPj. But to order P 4 , since each A brings 
at least one momentum power (the zeroth power is ac- 
counted for already in T), only combinations of the type 
(P1P2XP3P4), (P1P4HP2P3), (PP3XP2P4) can appear. 
The coefficients of the first two terms have to be equal be- 
cause of the cyclic symmetry of eqn. (|31H The coefficient 
of the last term is in general independent. 

In terms of fictitious momenta, all flowing into the di- 
agram, whose conservation law is Pi + P2 + P3 + P4 = 0, 
we obtain 



(40) 



The two numbers d\, d 2 contain the non-trivial infor- 
mation in this diagram. We have explicitly pulled out 
the color factor (3 as will be shown shortly) from the d's, 
which contain in this way only momentum and spin (since 
flavor will be dealt with at the end when the external legs 
are matched to the physical particles). 

The third term in l|4(JI) , to order 4 in momentum, with 
no ladder, is simply a wavefunction overlap given by the 
usual Feynman rules (notice an extra (-1) will be due to 
the fermion loop) 



J 



d 4 q 

(2tt) 4 \A 2 q 2 -B 2 



Tr 



A^(P u q)(A J q + B)A^(P 4 ,q)(A jq + B)A^(P 3 ,q)(A jq + B)A^(P 2 ,q)(A A + B)} (.41) 



The Dirac traces can easily be computed with FORM. 
The integral is then reduced by using tensor identities 
like 



F{q 2 )q tl q u q p q a = 



g^gP* 



9^9"" 1 4 



24 



to a one dimensional expression which can then be nu- 
merically evaluated once a specific model (and hence 
Bethe-Salpetcr wavefunctions) is chosen. Notice that this 
simple momentum routing is correct only to order P 4 . 

The first and second diagram in (|40|l contain a ladder. 
A simple way to calculate them is to write an integral 
equation for the object 



external pion momenta is 

U := U (fc 2 )P 1 -P 2 + U 1 (fc 2 )A;-P 1 fc-P 2 ^43) 

U 2 (fc 2 )P r P 2 / k + U ?> {k 2 )k-P 1 k-P 2 /: + U 4 (k 2 )k-P 2 ft + 

U 5 (fc 2 )fc-P! P 2 + U 6 (fc 2 )fc-P 2 jk ft + U 7 (k 2 )k-P 1 jk P 2 + 

Ug^ 2 ) P a / P 2 + U 9 (fc 2 ) Jk Pi /P 2 . 
The functions are obtained by projecting this linear 
integral equation (analogous to the Bethe-Salpeter equa- 
tion) 



iU (Pl,P2,fc) 



Ai 



(44) 



iU : = 



x- 
x- 



Ai 



(42) 



whose most general expansion up to second order in the 



with the matrix projectors I, /k,/k /Pi, ... /k /P\ /P 2 , 
which provides us with a linear system of eight integral 
equations for the Uj. Defining a convenient quantity D := 
2fc-Pi k-P 2 — k 2 P\-P 2 the projections are 



J 



4P(P r P 2 U +k-P t fc-PaUi) 
4P(P r P 2 U 2 +k-Pi fc-P 2 U 3 ) 
4Pi.P2Pfc.P2U4 



Tr[{2D - 2fc-P )t jP 2 + q 2 Pi P 2 )U] 

Tr[(-2P 1 -P 2 fi + 2k-P 2 Pi+ ft P! P 2 )U] 

Tr[(2k-P 2 Pi_.Fi ^-2(fc-P 2 ) 2 Pi +D P 2 - fc-P 2 jk ft P 2 )U] 



(45) 
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4P 1 -P 2 Dk-P 1 U 5 = Tr[(D ft - 2(k-P 1 f + k-P 1 ft ft ft)U] 

4P 1 -P 2 Dk-P 2 U 6 = Tr[{P 1 -P 2 ft ft-k-P 2 ft ft)U] 

4P r P 2 £>fc-PiU 7 = Tr[(-2k-Pi P V P 2 + P V P 2 ft ft + k-Pi ft ft ft)U] 

AP l -P i DU s = Tr[(k 2 P 1 -P 2 - k-P 2 ft ft + fe-Pi ft ft - k 2 ft ft)U] 

4iVfl,DUg = Tr[{Pi-P 2 ft-k-P 2 ft + k-P, ft- ft ft ft)U] 



and the inhomogeneous part of the equation can easily be 
written down from the first RHS diagram in l|44|) . Once 
the equation for U is solved in a computer, the diagram 
can be closed from the left to give 



A 4 



(46) 



and calculated as a simple integral. Here the color factor 
of 3 in H4(Jf) can be easily seen, since all three vertices are 
color singlets and carry 5 cc i in color space.. 



C. Contribution with T — 3 A. 

We will reduce the 3 — A contribution fixing the indices 
r 4 A 4 A 2 A 3 (the other permutations can at the end be 
easily generated). Employing (|25[l in 



r 4 AAi 



A,1 t A 



r 4 



Ai 



r 4 Ai 

- — * m — - 



^ A 2 



(47) 



(in the first diagram fix j = 3 to substitute l|25|l. in the 
second diagram employ j = 1) to obtain 





r 4 




A x 




— <-*< — 
— ►** — 








s- 1 




A 2 




s- 1 




A 2 




— <-»< — 
— >-•* — 








r 4 




A 3 



(48) 



r 4 Ax 

A 3 A 2 

Next one can apply (|22(l to the single T appearing in 
this expression to generate, after some simple manipula- 



tions: 



75 



v c 3 i 



Ai 

i 

A 2 
A 2 

A 3 



(49) 



From this expression, the two terms with Cj are zero in 
the chiral limit (as appropriate to this paper). The Cj's 
diverge for low energies in the chiral limit: they contain 
the pseudoscalar ladder pole, But in this diagrams, the 
object to the right of the ladder contains a product of 
2 A's, each of negative parity, the result carrying pos- 
itive parity. By the symmetry breaking pattern of the 
theory, no massless scalar, pseudovector or tensor meson 
pole can make the ladder divergent at low momentum. 
Therefore the terms with Cj can be discarded and we 
have to consider only (substituting ja — -775 , good in 
the chiral limit): 



Ai 



WW* 



A 2 



A 2 
A 3 



(50) 

Therefore our next problem is to evaluate diagrams such 
as: 



k Ai 



h-Pi-P* 



A 2 



(51) 



3d 3 Pi ■ P 3 Pi • P 2 + 3d 4 P 2 • P 3 P 4 • P 2 



which again explicitly displays the color factor and where 
the new constants d 3 and d 4 have to be calculated in a 
specific model. 

To reduce the ladder in this diagram we could attempt 
to use the Vector Ward Identity for the 7 M on the left, but 
the momenta flowing in the adjoining propagators would 
require this 7^ to be contracted with — Pi — P 2 = P 3 + P 4 
and not just with P 3 . Or in the right part of the diagram 
we could reuse our result for the two- A (U) vertex from 
the previous section. But again the momentum flow is 
not adequate. The solution to this impasse is to use both 
ideas, but in a momentum expansion. The ladder in this 
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diagram can be substituted by its momentum expansion 
(|26[l . Since there are three powers of momentum already 
committed (one is the explicit P 3 , the other two need 
to be one in each A) only one more power is needed. 
Therefore we can use the ladder expansion to order one, 
and diagram H51(l can be rewritten as 



0(P) 



0(P 3 ) 



(52) 



/flC 



k r-i k' 1— 1 



k 

0{P) 



(fc'+Pi+P 2 
O(P) 



ri 1 

(i)A 2 
0(P 2 ) 



Now we can use the vector Ward-Takahashi Identity, 
which is satisfied to 0(P^), on the first diagram and 
on the left ladder of the second diagram, allowing us to 
substitute 

/ P 3 ^iS-\q)-iS- 1 (q-P 3 ) . 

The matrix object with a ladder and two deltas, to 
O(PW), which appears on the second diagram, is just 
U as defined in (|42[l . Now it is straightforward to show 
that (|51|l is equal to 



(2kP 3 (A' jk-B')+Afo) 



( 



\ 



fc5 

(^fPi+p 2 ) 



A 2 



(q+-P 1 +P 2 ) (1) 



(53) 

where the diagrams have to be evaluated to order (P^) 
since an explicit power of P has already been used. The 
left diagram in particular gives rise to four simple sub- 
diagrams since two powers of P are committed in the 
A's, but the other power can be distributed alternatively 
between these A's or the two propagators which carry a 
power of P. 

These diagrams can all be evaluated easily as a sim- 
ple loop integral in the computer to obtain <i 3 , d 4 . The 
contribution from l|47|l is finally 

-(d 3 (Pi-P 3 P r P 2 - Pi-P 2 P 2 -P 3 ) 



( 

+d i (P 2 -P 3 P r P 2 - P r P 3 P 2 -P 3 ) + permutations) 



(54) 



D. Contribution with AITT. 



The contribution 



A x 



c 



r 4 



tr 2 



3 - C 



r 4 



Ai 
1 



r 2 



r 2 



(55) 

can be reduced as follows: 1) apply l|2U|) and (|23[) to the 
Ai in the first term, and (|12|l to the T 3 in the second 
diagram to obtain an expression similar to 148|) in which 
r 4 is isolated. 2) Employ l|22l) to eliminate r 4 . 3) Repeat 
the operation to eliminate T 2 . Obtain 



C&P1+P4 



<ffPl+P 2 



+permutations 



(56) 



In this expression, two explicit powers of P are present, 
namely P\ and P 3 in the vertices. The other two powers 
have to be produced from a propagator expansion. This 
is the third (and last) different diagram type that we need 
to parameterize: 



JL 



(57) 



3d 5 Pi • P 3 Pi • P 4 + 3d 6 P 3 • P 4 Pi • P 4 



Next we show how to calculate d§, d$. Since there is a 
ladder which contains powers of P, we need to recall the 
ladder expansion to second order in (|26|) . Eliminating 
the loose rung with the help of (|30|l . and employing the 
vector Ward Identity to generate a vertex 

V(P) = 2q.P(A'{q) & - B'{q)) + A(q) /P , 

we can show 



V(Pi) f~^~> ^(P3) + V-(Pi) f~^~rH~> V(P 3 ) - V(Pi) f ' Sfef^ V(Pa) 

J57j = ( ^ A+P ^ )<2) (^p 1 +p^P(mpZp^ (^p 1+ P4) (1) (^p 1+ P4) (1) (58) 

I 

The first and third diagrams are again straightforward a ladder. We can interpret this ladder as "dressing" ei- 
traces and integrals. Only the middle diagram contains 



9 



ther of the vertices, and write immediately an integral 
equation for one of them. Taking for example 

k 



(fe+pl+p 4 )' 



= V 



to the in the second, then using (|13f) to simplify the 
remaining T4 , neglecting the positive parity ladders when 
they are divided by a Cj containing the pion pole, re- 
absorbing the ladders and simplifying, one obtains: 



the function vertex V so defined satisfies a linear inho- 
mogeneous integral equation (the first argument is the 
momentum entering the diagram through the vertex, the 
second is the relative and the third the total momentum 
between the fermion lines at the vertex): 

V (p 1 ,fc + ^i-^,p 1 + p 4 ) = 

[2kPi(A'(k) fc-B'(k)) + A(fc)/Pi]S(fc+Pi+P 4 ) (1) S'(fc)- 1 



g 1— 1 , . ^1 



(64) 



/Pi 



+ per mut at ions 



V € 



q 



(60) 



As is evident, V admits an expansion up to second or- 
der in momentum identical to (|43[l in terms of a new 
set of functions Vo(fc 2 ), ■•■ ,Vg(fc 2 ). The integral system 
of equations i|6U|) is very similar to l|44|) . the only differ- 
ence being the inhomogeneous term. Therefore the linear 
projections in ll4"5|l still apply, and both systems can be 
solved with basically the same iterative computer code. 

Finally the middle diagram in l|58|l can be closed to 
read 



which can be written down immediately in terms of the 
d's defined in cqns. i|5T)l . (|57Jl as 

3[-d 3 P 1 -P 3 P 1 -P 2 - d i P 2 -P 3 P 1 -P 2 (65) 
+d b Pi ■P 2 P 1 -P 4 + d 6 P 2 ■ P 4 P 1 ■ P 4 + permutations] . 

f. ArAr 

Two permutations contribute: AiI^AsI^ and 
TiA2r3A4. The reduction is in all analogous to the pre- 
vious ones, yielding a contribution 



-3 



(V2U) 



(ffpl+p 4 ) ( 



■y V{P3) 



(61) 



- [d 5 {Pi -P 3 Pi-Pa + Pv P3 Ps-Pi + permutation) 

+d 6 (P 3 -Pi Pi- P4 + Pi- Pi P3-Pi + permutation)] 

(66) 



which is easy to evaluate with the help of a symbolic 
manipulation program. Finally we give the expression 
for Ij55(l in terms of the d's: 

3d 5 P r P 3 (P v P 4 + P r P 2 )+ (62) 
3d 6 {P 1 -P 2 P 2 -P 3 + P 1 -P i P 3 -Pi) + permutations . 



E. Contribution with AAIT. 

With two A corrections, there are two topologically 
distinct diagrams that can contribute. They are different 
because while reading around the fermion loop one can 
find the external legs in the orderrrAA or in the order 
TArA . We start by the first term, namely 



G. 4 — T contribution. 

The last piece stems from the term with four powers of 
r. By repeated use of l|22|) it can be shown to contribute 

-[d 5 (P 1 -P 3 P 3 -P i + P 1 -P 3 P 2 -P 3 ) 



+d 6 (P x -PiPs-Pi + P,- P 2 P 2 ■ P 3 )] . (67) 
(Here the cyclicity of l|31l) can be recovered by using 

E^ = o). 

Combining the results from subsections A through F 
and summing all permutations l|32|l with the isospin fac- 
tors H34I I35|l . the standard amplitude A(s,t, u), takes a 
form identical to (@J and we deduce 



tA 2 



Ax 

3. 



Ax 

V^A 2 



h = —{-2d 1 +d 2 -6d 3 + 2d 4 + 5d 5 +d 6 ) 
3 

h = - — (d 2 -2d 3 + 2di + d 5 -d 6 ). 
16 

(63) Therefore, to obtain the Vs numerically, one must: 



(68) 



(and the corresponding three more permutations) . Again 
by applying l|23|l to the Ai in the first diagram and i|l(J|) 



1. Solve the Schwinger-Dyson equations for the prop- 
agator, equation Q. 
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2. Employ the obtained A, B functions as input to 
the Bcthe-Salpeter equations (|16|) and solve them. 

3. Then use the obtained Fq, Go, ■■■ as input for equa- 
tions 1331) and (JHOJ to obtain U and V. 

4. Perform the integrals igBJl, C3J), (ED- 

5. Assemble gDJ, 0, @, <E3, ©, E3- 

In this external momentum expansion all integrals and 
integral equations are functions of internal variables only 
of quark momenta q 2 , k 2 , k-q. The diagrams could alter- 
natively be evaluated on the lattice. 

V. MODEL EVALUATIONS. 

We would like to provide simple model evaluations of 
all these calculations, well aware of model limitations and 
that a thorough phenomenological analysis can only be 
carried out with more sophisticated interactions such as 
employed in 0,0]- We employ two Feynman gauge mod- 
els featuring an interaction 

VK(q)V = ^K(q)^ . (69) 

This simple choice of a vector- vector interaction (as op- 
posed to the more popular Landau gauge transverse ten- 
sor kernel) simplifies the Gamma matrix traceology (in 
this calculation carried out with the help of two indepen- 
dent computer codes, one written in MATHEMATICA 
and one in FORM 22]) so that the standard Llewelyn- 
Smith BS wavefunction for the pion [7J reduces to 

X (P, k) = 75 (E(P, k) + F(P, k)/P + G{P, k) Jkk-P) 



= 75 ^E0(k 2 ) + y -^-E2(k 2 ) + F0(k 2 ) P 
+ GQ(k 2 ) jkk-P + ■ • •) (70) 

where the standard H function decouples from the rest 
of the system and therefore we further ignore, and the 
momentum expansion shown is complete up to second 
order for a symmetric momentum routing x (that is, 
the fermion lines out of x carry k + P/2 and k — P/2). 
The E0 term is determined by chiral symmetry to be 
r A {P-o,k) rp^g reg ^. Q £ k e gxpansion constructs the func- 

V2j/^ 

tion A(P, k). The same power series l|70|l can be written 
down for the axial vertex, 

T A (P, q ) = [(A(q-P/2)U-^)- (71) 
A(q + P/2){jq + f ) - (B(q - P/2) + B(q + P/2))]f 

expanding in powers of P, and up to a normalization we 
recover the equivalent to 170(1 



E0 A = 2B(q 2 ) (72) 
F0 A = ~A{q 2 ) 
G0 A = -2A'{q 2 ) 
E2 A = 2B"{q 2 ) 

Subtracting (|72|) from (|70|) we obtain some new functions 
of q which provide the needed expansions for A: 



A^(q + P 1 /2,P 1 ) = A^(q,P 1 ) = F0(q 2 ) f l +m{q 2 )q-P l h (73) 
M 2 \q + P 1 /2 1 P 1 ) = E 2 (q 2 ) { -^^ + FQ'(q 2 ) /Piq.Px + G0'(q 2 ) ^q-P,) 2 

I 



(valid for symmetric momentum routing when El, Fl, 
Gl all vanish). Further, with the 7^7^ kernel another 
trace-related simplification occurs in H44[) and (|60|) . and 
the functions U 6 , U 7 , U 8 , V6, V7, Vs equal the inhomo- 
geneous term in their respective equations, the homoge- 
neous (integral) parts of the equations being zero. 

All that remains is to consider some specific form for 
K(q 2 ). We will look at two models whose euclidean an- 
gular integrals can be done analytically, leaving only one- 
dimensional integral equations to solve numerically. 

The first is a simple Gaussian kernel (whose Euclidean 
angular integrals are Bessel functions J23j) 

K(q) = g 2 exp(-q 2 /A 2 ) (74) 



where g provides the coupling strength and A the scale 
of the interaction, (the results are sketched in table P). 
The second is a rational kernel 

K(q)=g 2 [l/(q 2 -\ 2 )-l/(q 2 -A 2 )} (75) 

where g, A and A represent some quark model param- 
eters fit to provide a good condensate and constituent 
quark mass via the Schwinger-Dyson equation. The eu- 
clidean angular integrals are straight-forward using rela- 
tions such as 

f 1 r^— 

/ = 7r(a — V a — 1) , 

J-i a-x 
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TABLE I: Results for the toy gaussian model in the chiral 



limit. 


Dimensionful mag 


nitudes are in 


MeV. 






A g 


M{k 2 = 0) 


-(**} 1/3 


u 


h 


h 


500 


6.3 


385 


222 


76 


-0.02 


0.066 


600 


6.3 


462 


267 


91 


-0.02 


0.060 


500 


6.5 


468 


237 


83 


-0.018 


0.062 


800 


5.5 


125 


213 


51 


-0.019 


0.080 



TABLE II: Results for the toy rational model in the chiral 
limit. Dimensionful magnitudes are in MeV. 



A 


A 


g 


M(k 2 = 


= 0) 




U 


h 


h 


250 


300 


13.0 


559 




242 


82 


-0.011 


0.11 


300 


350 


13.0 


471 




240 


77 


-0.015 


0.11 


500 


550 


14.2 


269 




247 


68 


-0.015 


0.12 


700 


770 


13.4 


155 




250 


56 


-0.012 


0.12 




L 


X-y/l — X 2 


n(a 2 


- 1/2 -a 




-1) ■ 








a 


— X 







Hence all integral equations are one-dimensional in mo- 
mentum space. 

When the dimensionful parameters, of order the strong 
interaction scale, are close in value, and for large enough 
g, then the potential is infrared enhanced supporting chi- 
ral symmetry breaking. Notice that the Rainbow-Ladder 
fermion loop diagrams constructed with these interac- 
tions are finite due to the exponential or l/q 4 high en- 
ergy behavior and no renormalization program is needed. 
This is associated to a gluon propagator scale, which de- 
termines the scale of the Bethe-Salpeter wavefunctions, 
which in turn regulate the meson loops that would appear 
in chiral perturbation theory. We are currently investi- 
gating this issues. 

VI. RESULTS AND DISCUSSION 

We have worked out the pion-pion scattering ampli- 
tude at low energy in the chiral limit to 0(P 4 ) from 
a general microscopic quark Schwinger-Dyson approach. 
Upon comparison with the chiral lagrangian formalism, 
the main result of this paper is a pair of relations that 
would in principle allow to directly evaluate the coeffi- 
cients 1%, I2 in any specific model (i.e. after choosing 
the Lorentz structure of the quark-antiquark interaction, 
strength and shape of any potential or dressed gluon 
propagator) provided it supports the standard mecha- 
nism of chiral symmetry breaking |l7j , |l8j , p3 | . symme- 
try properties of the Y pion vertex, are model dependent 
and are of order P 4 (therefore vanishing at low energy) . 
Lorentz invariance restricts their form and allows only 
for two coefficients. 

On a first glance, diagrams (|4(JI etc.) seem as difficult 



to calculate as the full pion-pion scattering amplitude, 
but one needs to remember that each A pion wavefunc- 
tions entering the calculation can be taken to have only 
one power of the external P to this order, since there 
are four of them and they vanish by definition at P = 0, 
and the propagators can also be taken at P = 0, leaving 
just one momentum k around the loop. This is a major 
simplification. 

In the previous section we have shown how the numer- 
ical evaluation can be performed with two very simple 
(too simple) kernels. Our results for li and I2 have phe- 
nomenologically correct signs and ratios, specially for the 
gaussian toy model, see tables HI and ITT1. Possibly because 
our results are computed at P = and without pion 
loops, our l\ and I2 seem too large in absolute value. This 
can be appreciated SIEISIlSililiilllllliilP in 
tables InTl II VI and Ivl which summarize the present status 
of knowledge of these coefficients. The various conven- 
tions used in the literature (SU(3), SU(2), renormalized, 
barred, etc.) have been unified to the renormalized I's at 
the scale of the p meson. Some of those determinations 
carry information relative to finite quark (and therefore 
pion) masses, to kaons and etas, or to pion loops, all of 
which have not been taken into account in this work. The 
fairest comparison therefore is with the oldest results of 
Andrianov ^(J in the large N c limit (no gluon correc- 
tions) and Pham and Truong [27). The later authors 
obtain the interesting relations 

f 2 
m z p 

and 

1 f 2 
6 m% 

If the a mass is sent to infinity, then the large N c ratio is 
recovered. For finite a masses between 350 and 750 MeV 
we obtain the range given in tablc lTvl This demonstrates 
the importance of repeating the model evaluations with 
kernels whose meson excitations in various channels are 
known. 

Our approach is potentially superior to the resonance 
saturation approximations since it includes the full ver- 
tex and ladder structures, that is, effects of continua and 
higher resonances. In particular we include the four pion 
direct interaction, and the exchange of the full series of 
excited vector and scalar mesons. We also stress that the 
masses of our p and a mesons are expected to be of the 
right order of magnitude, because our constituent quark 
mass have reasonable values of the order of 300M eV to 
AQQMeV, see tables|I]and[n] Our approach is also poten- 
tially superior to Nambu- Jona-Lasinio determinations by 
allowing to lift the approximation of contact interactions 
between quarks. Thus we pave the way to calculating the 
li coefficients from the lattice or from accurate Schwinger- 
Dyson solutions. Finite current quark masses and meson 
loops remain to be incorporated to improve the precision 
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of our calculation. We are currently considering some of 
these issues. 



TABLE III: Phenomenological determinations of the 1 pa- 
rameters (fits to scattering lengths or phase shifts in pion 
scattering.) 



Authors 


Z[(m p )-10 3 


l r 2 (m p )-10 3 


Gasser & Leutwyler 


-4.2 ±3.9 


9.0 ±2.7 


Bijnens, Colangelo, Talavera; 






Colangelo, Gasser, Leutwyler 


-2.2 ±0.6 


9.0 ±2.7 


Yndurain 


-4.1 ±0.7 


9.8 ±0.5 


Gomez Nicola & Pelaez 


-3.3 ±0.7 


4.8 ±0.6 



TABLE V: Theoretical determinations of the 1 parameters 
(based on the large N c approximation and/or the Nambu- 
Jona-Lasinio model). 



Authors 


Z[(m p )-10 3 


Z£(m p )-10 3 


Espriu, de Rafael, Taron 


-3.18 


6.3 


Bijnens, Bruno, de Rafael 


-4.8 


6.4 


Ruiz Arriola 


-7.4 


7.8 


Andrianov 


-3.2 


6.4 



TABLE IV: Phenomenological determinations of the 1 pa- 
rameters (Based on p meson resonance saturation). 



Authors 


Zi(m p )-10 3 


^(m p )-10 3 


Gasser & Leutwyler 


-8.4 


8.4 


Pham & Truong 


-5.5 to -24 


14 


Ecker et al. 


-6.1 ±3.9 


5.3 ±2.7 


Pennington & Portolcs 


-2.1 ± 1.2 


5.9 ± 1.0 
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f^-y d^S 1 -T A £^y d » s o 1 75-p 



(A5) 



re-absorbing the remaining ladder, and employing the 
free propagator and bare axial vertex in eq. (0 this is 
proportional (because of the isospin factor not included) 
to 



Tr J (i ^sd li s- 1 s-r A _ I ,s 



7 M 75 
2i 



S) 



(A6) 



APPENDIX A: NORM AND /„ 

In this appendix we present the proof of the well known 
fact that the normalization of the Bethe-Salpeter wave- 
function is the pion decay constant, rewritten in terms of 
our Ward Identity techniques. At null external pion mo- 
mentum, \ a (P = 0) and V\ are proportional because eq. 
(fTUfl reduces to the homogeneous BS eq. l(TB|) . Introduce 
an arbitrary proportionality constant n x by means of 



X a (P = 0) 



pa 



(Al) 



the i guarantees \ t° be real because of the explicit form 
of T in eq. 

The normalizing condition for the Bethe-Salpeter so- 
lution is, following Llewelyn-Smith 30], 



2iP fl 6 



ab 



(A2) 

Because the right hand side of the normalization condi- 
tion is of first order in P M we can expand in X and in Ta 
as in (|2(J[1 . There are two terms, depending on whether 
the derivative is applied to the upper or lower propaga- 
tor. 

The first term can be written as 



r.4 



T A a _, 



X-p 1 




(A3) 



The term with IT can be shown to be zero. To see it, 
one needs to take a derivative of eq. J7J which gives 



sd^s = ^Q^y^o 1 



(A4) 

which applied to the FT term in l|A3[) and employing eq. 
reduces it to 



Eliminating T with its Ward Identity l|13H , and after some 
elementary operations, this equals 



iTr / (j u d^S - d v S^)P u = 



(A7) 



Returning to equation i|A3(l and once the term with IT 
has been shown to vanish, we need to evaluate the terms 
with a x and a T. 

Diagrammatically again, one has (recall equations 
(El, (E3 and JUJ) 



1 A_ 



fc+P/2 

--Tr{ l5 _ p ^ Xbp d^(k + P/2)} 
~2i 



(A8) 



We get two terms. The first is of order P 2 because of 
the c in the denominator, and vanishes. The second is 
indeed non zero. Going back now to eq. (|A2|) . we obtain 



2iP„ = 



.757/.— 



(A9) 



The definition of the weak decay constant, /„• yields (with 
no pion loops) 



X b <^>757m^" 



(A10) 



and therefore we must have 



-IV i 



-71 



p 2 



Ta" 



W"ir fir 

which immediately leads to equation (|2U[1 . Finally, direct 
calculation of equation IjAlOfl leads to 
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if, = 3 



d 4 q 



(2tt) 4 (A 2 q 2 - B 2 ) 2 



EO 4AB + 2q B 



dA 

dW) 



-A 



dB 
W) 



F0{2q 2 A 2 + 4B' 2 ) + q z GO{B l - q z A 



2 a2\ 



iAU} 



Notice the explicit color factor of 3. All through the pa- 
per the Bethe-Salpeter wavefunctions have been taken 
proportional to the identity <5 CC < in color space. Had 
we normalized them in a different way, say 5 cc i /V3, this 
would immediately affect the formula above reducing the 
factor to \/3, the rest being absorbed by the functions EO, 
FO, GO. The low energy theorems (Gell-Mann-Oakes- 
Renner theorem, Weinberg's amplitude, etc.) are un- 
changed by this choice since the explicit form of these 
functions is never used to prove them: they are always 
eliminated in terms of f K . But the l\, I2 constants of 
the chiral lagrangian would indeed have to be rewritten 
in terms of the modified Bethe-Salpeter wavefunctions. 
This of course would not affect its numerical value. 
The normalizing condition (|A2|) . 

X-p<" ( ) >Xp 

^ a P M V -~ J J = 2iP fl . (A12) 

which yields (|A10|I can also be directly evaluated without 
using Ward Identities. By taking a derivative of this 
equation respect to P M (and contracting over /j, as usual), 
we get the following (derivatives respect to P act only 
on the function immediately behind them and the color 



factor is explicit). 

8 ^ 3 /(2^ Tr[ (M3) 
2d„xAP)9^S(q + P/2)xA-P)S{q - P/2) + 
2 Xv (P)d^S(q + P/2)xA~P)S(q - P/2) + 
2 x ,(P)d^S(q + P/2)d»xA-P)S(q - P/2) + 
2 Xv (PWS(q + P/2)xA-P)d„S(q - P/2)+] 



To calculate this normalization one needs to make the 
wavefunction explicit 

XAP) = To(E0(q 2 ) + F0(q 2 ) P + G0(q 2 )q-p k + • • •) 

(A14) 

which evaluated at P = yields 

XAO) = l 5 E0(q 2 ) ; d^xAP = 0) = j 5 (F0 ltl +G0 h%) . 

and make use of the propagator expansion, (|27l 1281 129|) 

above. A simple check on the resulting expression is to 
substitute x by T, which yields 



d A g 
(2tt)4 



-2B(d p B)d p 



A 2 q 2 -B 2 ) B dpd " ^A 2 q 2 - B 2 



= 



which vanishes upon employing Green's first identity. This checks the zero in (|A7|) with a completely independent 
calculation, and is also approximately observed in our computer codes. The result for is 



3i 



d 4 q 



1 



(2tt) 4 A 2 q 2 - B 2 



[(AB' -BA , )2q 2 B(Fo + q 2 G ) 
-AB 2 (4P + q 2 G ] 



(A15) 
(A16) 



which, upon comparison with. (|A11() provides an integral constraint between the Schwinger-Dyson and Bethe-Salpeter 
solutions. The barred quantities have been defined in i|73|) . 

I 



To conclude this discussion we recall the derivation 
of the Gell-Mann-Oakes-Renner relation in the Bethe- 
Salpeter formalism, this has been presented in 0], [l4| . 
together with the discussion of the Weinberg theorem. At 
zero external pion momentum the Axial Ward Identity 
reads 



and as previously discussed 



(Ais) 



the vertex to write 



2imTl{k,k) 



(A17) 



To start the simple demonstration of GMOR, "undress" 
te 



ST%S 



lb 2 



(A19) 
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and, neglecting the contribution of higher pion states 
(which exactly decouple in the chiral limit, as in ref. |24[) 
we can saturate the ladder by the pion pole yielding 



substituting now eq. (|A1|) in the form 

a _ s^^ + r's- 1 <J a 

X - o 

and comparing with eq. I|A18() we immediately obtain 



(A20) 



(A21) 



corresponding to the GMOR relation upon identification 
of = f w , consistent with the Llewelyn-Smith normal- 
ization condition. 

We finally remind the reader of the explicit expression 
for the BCS planar condensate: 
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A k jk- 



= TrS 
Bk 



A k k*- 
d 4 k f 



Bl 



+ is 
4B k 



(2tt) 4 A\tf 
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(A22) 



-2mTrS = n\Ml 



